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Abstract 



We consider an electron-phonon system in two and three dimensions on square, hexagonal 

and cubic lattices. The model is a modification of the standard Holstein model where the 

'^ ■ optical branch is appropriately curved in order to have a refiection positive Hamiltonian. 

Q . Using infrared bounds together with a recent result on the coexistence of long-range 

CJ \ order for electron and phonon fields, we prove that, at sufficiently low temperatures 

and sufficiently strong electron-phonon coupling, there is a Peierls instability towards a 

period two charge-density wave at half-filling. Our results take into account the quantum 

fluctuations of the elastic field in a rigorous way and are therefore independent of any 

O ! adiabatic approximation. The strong coupling and low temperature regime found here 



is independent of the strength of the quantum fiuctuations of the elastic field. 
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1. INTRODUCTION 

The Holstein model was originally introduced [Q to describe the motion of polarons 
in molecular crystals. Despite the simplifications introduced, it contains the essential 
features necessary to describe the interaction between lattice vibrations and itinerant 
electrons in a crystal. In the simplest version, the elastic field is modeled by a collection 
of Einstein oscillators with a single frequency associated to each site of the lattice and 
the electrons are treated as noninteracting fermions. This is a situation where itinerant 
electrons interact with a fiat optical phonon branch. More sophisticated versions of 
molecular crystal models take into account a dispersion in the optical branch, as well 
as the Coulomb interaction between electrons. When a single electron is present (or 
when the electronic density is low), the Holstein model has been widely used to study 
the composite entity formed by the motion of the electron and its associated lattice 
deformation, the so-called "Holstein polaron". We refer to 0, ^, ^ for recent works on 
the subject. 

Another situation of interest is the one where the electron density is large and where 
instabilities such as superconductivity (SC), charge density wave (CDW), spin density 
wave (SDW) occur. In the present work we focus on CDW formation for the Holstein 
model at half-filling. The prediction that such an instability may be present in electron- 
phonon systems goes back to Peierls and Frohlich 0. They showed within the 
adiabatic, or more accurately static, approximation for the elastic field that, for one 
dimensional electron-phonon systems, a periodic modulation of the lattice commensurate 
to the electron filling would open a gap at the Fermi level. Remarkably, the electronic 
energy lowering caused by the gap opening is larger than the elastic energy cost and 
therefore a "Peierls instability" occurs. This phenomenon plays an important role in 
organic compounds such as polymer chains ||^ and is for instance responsible for bond 
alternation in polyacetylene chains (CH)^; |Q] or annulenes. In two or three dimensional 
systems, this instability is not generic because of the more complicated shape of the 
Fermi surface. However, if the Fermi surface has nesting properties (for example on a 
square lattice at half-filling it is perfectly nested), the Peierls instability will occur for 
the same reasons than in one dimension. 

The analysis of Peierls and Frohlich and most of the subsequent work has been based 
on the adiabatic approximation for the elastic field. In many respects, it is still an 
open problem to determine when the Peierls instability survives against the quantum 
fiuctuations of the elastic field. According to the analysis of Hirsch and Fradkin and 
to perturbative renormalization group calculations |]10|, the Peierls instability should not 
occur for weak electron-phonon coupling if the itinerant fermions are spinless but should 
occur for any coupling when spin is taken into account (at least at half-filling). However, 



recently this view has been challenged by numerical work |jTT| which predicts that the 
instability is destroyed at weak electron-phonon coupling even when spin is involved. 
For strong electron-phonon coupling, it is likely that the CDW forms in all cases. 

There exist a few rigorous studies of the question starting with a result of Kennedy 



and Lieb [12| for the (static) Su-Schrieffer-Heeger (SSH) model. Their analysis shows 
rigorously that a period two ground state forms at half-filling consisting of alternating 
short and long bonds when electron-electron interaction is not taken into account. This 
result has been extended when a Hubbard interaction is added to the SSH Hamiltonian, 
by using a Jordan- Wigner mapping of the problem on a spin model [|T3|. It is shown 
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that at half-filling the ground state is either translation invariant or period two. The 
static Holstein model has also been the object of rigorous work in one, two and three 
dimensions. It is known that for the half-filled situation the ground state is period two 
in any dimension and that this symmetry breaking persists for low temperatures for 
dimensions greater or equal to two |T^. Away from half- filling, it has been proved in 
the one dimensional static Holstein model that the Peierls instability also occurs for 



rational densities |]T6|. When quantum fluctuations of the elastic field are taken into 
account Benfatto, Gallavotti and Lebowitz |T^ have shown, using a rigorous version of 
the renormalization group analysis valid for spinless electrons in one dimension, that 
there is no instability for weak electron-phonon coupling and the electron correlation 
functions are those of a Luttinger liquid. We note that their analysis holds for any 
electron filling but breaks down when electron spin is included. Also it does not cover 
the case of large electron-phonon coupling. A very general result of Freericks and Lieb 



15| asserts that for any even number of electrons on a finite lattice, the ground state of 
the Holstein model is unique and has zero total spin. This property is compatible with 
the presence or absence of a Peierls or other instabilities and with the uniqueness or 
non-uniqueness of ground or Gibbs states in thermodynamic limit, as long as the system 
is in a singlet state. 

In this work, we consider the two and three dimensional situation for a modified 
Holstein model taking into account quantum fluctuations of the elastic field exactly. 
Our model consists of non-interacting spinless electrons on square, hexagonal or cubic 
lattices interacting with an optical phonon branch which is appropriately curved. For 
this special phonon branch, we are able to use infrared bounds together with a recent 
result on the coexistence of long-range order for electron and phonon fields, to prove 
that a period two CDW forms at strong coupling and half-filling. Our result shows 
that at least for some regime of parameters the Peierls instability survives quantum and 
thermal fluctuations. The regime that we identify is independent of the strength of the 
quantum fluctuations. Our use of infrared bounds is based on the recent observation 
that some fermion Hamiltonians have the reflection positivity property if the hopping 



matrix satisfies some sign conditions [ff^, |T3 . These conditions are most easily expressed 



in terms of the "magnetic flux" through a plaquette, i.e. the phase of the product 
of hopping amplitudes around an elementary loop of the lattice. For a square or cubic 
lattice it turns out that the magnetic flux through squares must be equal to tt, whereas for 
hexagonal lattice the magnetic flux through hexagons must vanish. We note that it is not 
necessary to interpret this sign condition in terms of magnetic flux. For example in some 
situations the arrangement of orbitals of neighboring atoms may lead to overlap integrals 
having these sign properties. Such a situation has been encountered in coupled polymer 
chains |2^ and has also been the subject of speculations in the chemistry literature [pl[ 



The paper is organized as follows. In section 2, we give the formal deflnition of the 
model, the main results and we outline the general strategy of the proof. Section 3 deals 
with the coexistence of long-range order for electron and phonon flelds in a form which 
relates their associated correlation functions. This is then combined with the infrared 
bounds derived in section 4 to complete the proof for the square and cubic lattices in 
section 5. The necessary modiflcations for the hexagonal lattice are given in section 6. 
Finally we conclude by a discussion of a few open problems, that may be studied within 
the present formalism. 
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2. Definition of the model and main results 

We consider a system of electrons and phonons interacting on a finite hypercubic 
lattice N X . . .X N with N even, in d dimensions and with periodic boundary conditions. 
We call |A| = A^*^ the volume of the system. The Hamiltonian under consideration has 
the form 

H!, = Hl + Hf + H';^K (2.1) 

if^ is the kinetic Hamiltonian of spinless electrons 

Hl=Y. ^-y'^tcy (2.2) 

x,y£A 

where c^ and c^ are the creation and annihilation operators for an electron at site 
X. They satisfy the canonical anticommutation relations {c^,Cy} = {Cx,Cy} = and 
{c^,Cy} = 6xy The hopping amplitudes t^y = for |x — ?/| 7^ 1 and t^y = e^^^'y for 
|x — yl = 1 where the phases O^y are chosen in such a way that there is a magnetic fiux 
TT across each elementary plaquette of the lattice, i.e. ^piaquette ^^y ~ ^- ^^ choose A^ 
to be even and, as explained in the introduction, impose such an external magnetic fiux 
in order to have a Hamiltonian which exhibits the refiection positivity property. 
The Hamiltonian for the phonons H\ is 

Hf = Y,^{k){hth + 1/2) (2.3) 

k 

where the creation and annihilation operators for a phonon with wave vector k = 
(fci, . . . ,kd) satisfy commutation relations [b'l,b'l,] = [bk,bk'] = 0. All fc-sums are over 
the first Brillouin zone associated with the lattice. In the standard Holstein model one 
has a fiat dispersion relation uj{k) = cu = \J KJm where m is a mass and n the oscillator 
constant |^ . Here we consider a slightly curved dispersion relation of the form 



uj{k) 



\ 






+ — (rf + VcosA;,) (2.4) 






where a is some positive parameter, the reason being that ( |2.4| ) is compatible with 
refiection positivity. In fact there is a general class of dispersion relations that we could 
allow, the simplest of which is ( |2.4| ), but we will not pursue this issue here. 
For the electron-phonon interaction we choose the simple form 

Bt = 9Y.^^=={b^k-.bl)p^ (2.5) 

where p^ is the Fourier transform of the electronic density n^ — 1/2 = c^Cx — 1/2 
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It turns out to be more convenient to use the direct space language. We introduce the 
conjugate variables Sx and px with [si;,^^] = i5xy, 

-^= , (5fc + 6+1.) 

„ik-x 1 

Px = Yl -^=7^V^^Mk){h - htk). (2.7) 



lA|2i 
In these variables, the Hamiltonian for the phonons has the form 



^2 



< = E|| + fE^^+" E (^.+-.)' (2-8) 

x&h x&K \x-y\=l 

and the electron-phonon interaction is on-site 

Hr = 9Y.'x{nx-\). (2.9) 



xgA 



The Hamiltonian (|2.1| ) is invariant under the following unitary transformation Sx — > —Sx, 
c+ -^ (— l)'^'ca;, Cx —^ (— l)''*'c+. As a consequence, we have that the thermal expectation 
value (denoted by (— )a) of the number of electrons is half the number of sites and the 
one of the Sx variables is zero: {nx)^ = |, {sx)a = for any temperature P~^. 

In this work, we are interested in the two long-range orders for the Sx and Ux variables 
at ko = (tt, tt, vr) in three dimensions and /cq = (tt, tt) in two dimensions, i.e. we will show 
that there exists ei and €2 independent of A such that 

±-{s^,^Sk,)A = r^ E (-1)'''^'''(^x%)a > ei > (2.10) 

' ' ' ' x,yeA 

and 

T^Ap-koPko)A = r^ E (-1)'^'^'''((^^ - V2)(r^, - 1/2))a ^ e^ > (2.11) 

' ' ' ' x,yeA 

for sufficiently large A and some range of parameters. More precisely we show that long- 
range order is present on the square lattice at zero temperature if {K, + 32a)/g^/a ^ 0.19 
and at sufficiently low temperatures if (k + 48a)/g^/a ^ 0.22 for a cubic lattice. It is 
of interest to notice that the strong coupling regime obtained here is independent of the 
mass m of the oscillators. 

On the hexagonal lattice (section 6) we obtain a similar result for zero magnetic flux 
(i.e. all hopping terms have the same sign), namely ( p.lO| ) and (p.ll|) are valid at zero 
temperature for (k + 24a)/5'-\/a ^ 0.18. 

It has to be noticed that in the two dimensional cases (square and hexagonal lattices) 
our result is limited to zero temperature but this is probably an artifact of our technique. 
Since the symmetry breaking involved in this model is discrete (Ising like), we expect 
the result to be valid also at low temperatures. 

One would think that the adjunction of the a-term in the phonon Hamiltonian is able 
to create the long-range order (p.lO|) even at g = since it has tendency to force two 



nearest neighbors Sx and Sy to have opposite signs or, in the phonon language, the mode 
ko is the most favorable energetically since uj{k) is minimum for k = k^. This is in fact 
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not the case. Indeed we can compute explicitly in the case g = the expectation value 
we are interested in. We find 

showing that there is no long-range order since (|2.12|) tends to zero as |A| -^ oo. We 
can thus conclude that the long-range orders will appear only if the electron-phonon 
interaction is present. For this reason, this modified model provides an example of 
long-range order induced by the electron-phonon interaction in the presence of realistic 
quantum fiuctuations. 

Let us summarize the main steps for the strategy of the proof. First, in section 3, we 
establish the coexistence of the two long-range orders ( p.lO| ) and (|2.11| ). For this we use 



results developed in a previous paper |^ with which we can get a relation between the 
two-point Duhamel functions (s_fc, Sk)A and {p-k,Pk)A for any k. If we use this relation 
ai k = ko, this shows the coexistence we are looking for. However, in the next steps 
of the proof, we will also need this relation for k ^ k^. In section 4, we find an upper 
infrared bound for the Duhamel two-point functions (s_fc, Skjk using refiection positivity 
and, with the help of the previously mentioned relation, we deduce an infrared bound 
for {p-k,Pk)A- In section 5, we use a theorem of Dyson-Lieb-Simon |^ to obtain an 



upper bound on the two-point correlation function {p-kPk)A and then the sum rule 

|A 
1 



J2(P-kPk)A = E((^- - V2)^)a = ^ (2.13) 

k xeA 



to prove the existence of long-range order (|2.11|) for the electronic variables at fco. A 
more direct way would have been to combine the infrared bound on (s_fc,Sfc)A with a 
sum rule for the s^ variables: 



^x 



J2{s-kSk)A = J2^sI)a. (2.14) 

k xeA 

It would be sufficient to have a lower bound for this quantity, but we were unable to find 
a convenient one. We will come back to this point at the end of section 5. 

3. Coexistence of long-range orders for electron and phonon fields 

In a previous paper, we have developed a general formalism to discuss the possible 
coexistence of two long-range orders in a quantum system at finite temperature p2| 



which we summarize below. The Duhamel two-point function of two operators F and G 
is defined by 

(F,G)a = — f dvTi (e-^^^'^Fe-^^-'^^^^^G). (3.1) 

Za Jo 

The Duhamel two-point function {F~^,F)a is related to the symmetrized two-point cor- 
relation function through upper and lower bounds 



i(F+F + FF+)A/(/iA(F)) ^ iF+,F)A ^ ^(F+F + FF+)a (3.2) 
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where 



{[F\[pH^,F]])j, 

and the function f{u) is defined imphcitly for m > by 

/(■utanh-u) = — tanhw. (3.4) 

u 

The function f{u) is continuous, convex and strictly decreasing with hni„_»o f{u) = 1 

and hm^^oo f{u) =0. 

Let us suppose that the following commutation relation is satisfied between three local 

operators A^, B^ and Cx'- 

[HA,Ax]=fxBx + uCx (3.5) 

where n and u are two complex numbers. Let Ak, Bk, Ck the Fourier transforms defined 
as in ( |2.(:j ) . As a consequence of ( ^.5|) w e have the following relation between the Duhamel 
two-point functions of B^ and Ck |p2|| 

[\^i\./{Bt^^-\v\./(^^^^ ^ ^{[Al,Kk])K (3.6) 

where Kk = [HA,Ak]. 

With the help of inequalities (|3.2| ) and (|3.6| ) , we deduce the following general theorem 
of coexistence [^ : 



Theorem 1. Assume there exist three local observables A^, B^, Cx satisfying ( y7 
Suppose also that, for a given k, there exist three positive constants ak, bk and Ck in- 
dependent of \A\ such that ([y4+, [i^A, ^fc]])A ^ «fe; ([-B^, [-f^A, -Bfc]])A ^ bk and 
{[C^, [Hx,Ck]])A ^ Ck- Then for any temperature (3~^ and any sufficiently large A, 
there exists ei independent of A such that 

-^{B+Bk + BkB+)A > e,>0 (3.7) 

if and only if there exists €2 independent of A such that 

^AC^Ck + CkC+)A > 62 > 0. (3.8) 

To apply these results to our Hamiltonian ( p.l|) , we note the following commutation 
relation 

[H,Px] = i[{K + 8ad)sx + 4:a ^ Sy] + ig{nx - 1/2) = iSx + ig{nx - 1/2). (3.9) 

y 

\y-x\ = l 

Since the Fourier transform of Sx is 

d 

Sk= [k + Aa{2d + 2 J^ cos h)] Sk (3.10) 

l6D becomes 

y^/iS_k, Sk)A - gViP~k, Pk)A] < -^{\P-k,[H,Pk]])A- (3.11) 
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The double commutator on the right-hand side is simply a constant equal to K + 4:a{2d + 
2J2i'^oski). Finally, ( p.ll|) reads 



2 



i=l 



[k + Aa{2d + 2^ cos A;,)] ^/{s-k,Sk)A - g\/{p^k,Pk)j 

^ -[K + Aa{2d+2^coski)]. (3.12) 



P' 



i=l 



We will need inequality ( ^.121 ) in the next sections for k ^ ko. For the case under 



consideration in this section, i.e. to prove the coexistence of the two long-range orders 
at fc = fco, we evaluate (|3.12|) at A; = A;o to get 

I^V{S-ko,Sko)A ^ 9V{P-ko,Pko)A- J^- (3.13) 

We can now get information about the two-points correlation functions using the lower 
and upper bounds ( ^.21 ). Since the Sx are commuting variables 

{S-ko,Sko)A ^ {S-koSko)A- (3-14) 

For the lower bound involving the electronic densities, we need to compute the double 
commutator 

[n_k, [H, Uk]] = j-TT 5Z ^'■^y'^'^ ^°^[^ ■ ^y ~ ^)] ~ '^)'^tcy (3.15) 

' ' x,yeA 



The expectation value of ( |3.15| ) is bounded by 8d since \txy{2 cos[k{y — x)] — 2)| ^ 4 and 



{CxCy)\\ ^ 1. Since the function / is decreasing, we can replace the double commutator 
by 8d and get 

{p-ko,Pko)A > {p-koPk„)Af {-, ;-) (3.16) 

\{p-koPko)Aj 



where we also used the fact that the Ux variables commute. Finally, from (|3.13| ), ( p.l4| ) 
and (|3A6| ) 



Ky/{s_koSko)A > 9\ {p-koPko)Af 



\P~koPko)A 



This last inequality shows that if {p^koPko)A = 0(|A|), we have also that {s^koSko)A = 
0(|A|) since the function f{x) goes to 1 as a; goes to zero. In section 4, we will prove 
the existence of the long-range order for the electronic densities at fc = /cq- Because of 
( p.lTp , this implies the existence of the long-range order in the phonon variables. 

To close this section, let us notice that the commutation relation (|3.9|) in the case 
g = immediately leads to the absence of long-range order in the phonon variables for 
any k. This is consistent with the exact computation ( p.l2| ). 
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4. Infrared bounds 

The first step is to derive an infrared bound for tfie Duliamel two-point function of 
the phonon variables (s_fc, Sk)A for any k j^ ko hj using reflection positivity. We couple 
the phonon field in the Hamiltonian to a real symmetric field hxy, h^y = hy^, as follows 

2 

^A({M) = E£ + tE^'+" H {s. + Sy-h.,yf 

xgA xgA \x-y\=l 

+ ^ t,,4c, + f/^s,K.-l/2). (4.1) 

\x-y\=l xGA 

In the form ([4.1|) Hj^{{hxy}) is not reflection positive. However it can be transformed into 
a reflection positive form after a succession of appropriate transformations. We stress 
that these transformations cannot be done in general for fermions systems because the 
kinetic term usually does not have the correct sign properties, except if a flux vr is 
imposed through each plaquette of the square or cubic lattice. This has been shown 
in detail recently |]19[ and using these results one deduces in a standard way that the 
partition function Z^ associated to ( [4.1| ) satisfies Zf^{{hxy}) ^ ^a- 

We can now expand the left hand side of this last inequahty to second order in {h^y} 
to obtain 

( Y^ f^xySx, Yl ^^'?/'^^')a ^ ^ Yl ^^y (^-^^ 

\x-y\ = l \x'-y'\ = l \x-y\ = l 

This inequality has been derived for real symmetric h^-y but it is possible to extend it to 
complex numbers: 

( Yl ^*^y^^' Yl ^^V^xOa < ^ Y i^^s/i^- '^^•^^ 

|a;— j/|=l \x'—y'\=l \x—y\=l 

For each k in the first Brillouin zone, we choose h^y to be 

g-ifc-x _j_ Q-iky 
h^y = y= . (4.4) 



|A 
This immediately leads to the infrared bound for the phonon variables 

(s_fc,Sfc)A ^ — . (4.5) 

Apa [2d + 2 2^. cos ki) 

From the relation (|3.12|) we transfer the information contained in (|4.5| ) onto the 
Duhamel two-point function for the electronic densities 

9\/{p-k,pk)k ^ [k + 4a(2d + 2 ^ cos ki)] ^/{s_k, Sk)A 

i 

+ -[fi; + 4a(2rf+2j]cosA;,)] (4.6) 
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and inserting ( |4.5| ) 

r- — K + 4a(2rf + 2^,.cosA;,) /k + 4a(2rf + 2 ^. cos fc^^^ ^^^ 



< 



< 



K + 16ad 



P 



^A(3a{2d + 2^,cosh 



y/n + 16ad 

7^ 



K 



+ IQad + v^Kl^T6c^A/4a(2rf~+~2^~cos^ 



< 



2{k + I6ad) 



^/4p^42dT2Y~mski) ^4:(3a{2d + 2^^cosk 

In summary, we have also an infrared bound for the electron density 



(4.7) 



{P-k,Pk)A ^ 



/5a^2(2rf + 2^.cosfci) 
5. Proof of long-range order for the square and cubic lattices 



(4J 



Inverting the lower bound in ( |3.2| ) we obtain an estimate for the two-point correlation 
function of the electronic variables 



l/3'Ck 



{p-kpk)A ^ - \/BkCk coth J —^ 



(5.1) 



where Bu is defined in 



and Cfc is the upper bound of the double commutator ( ^.15| ) 
{[n_k,[H,nj,]])A ^ 8d = C^. (5.2) 

Using ( |5.1| ) and the sum rule ( |2.13| ) in the limit of large volumes, we get 

1 



1 , , 1 

rn-VP-fcoPfco/A ^ - 



coth 



^ V2{k + IQad) 



(2vr) 
-^^ lda{2d + 2Y,co.h 



d 
a{2d + 2 J2i cos ki) 



(5.3) 



Inequality ( |5.3| ) shows that there is long-range order whenever the right hand side is 
strictly positive. Because of theorem |I]the long-range order occurs both in ■ni{p-koPko)A 

and in rl-^{s_koSko)A- 

The integral in ( |5.3| ) can be computed numerically (when (3 = oo) which leads to the 
condition stated in section 2 for three dimensions 

^""^^"^ ^ 0.22, rf = 3, and large f3. (5.4) 

gy/a 

In two dimensions, ( |5.3| ) is not a convergent integral. One must first take the limit 
(3 ^ oo which yields an inequality similar to ( |5.3| ) but with the coth replaced by 1. 
Computing the remaining integral leads to the condition 

(k + 32a) 



gy/a 



^ 0.19, d = 2, and /3 = oo. 



(5.5) 



If we fix a and n, we see that there will be long-range order for large values of g. On the 
other hand, if we fix g and a, there will be long-range order for sufficiently small k. These 
two results are intuitively plausible: there is long-range order for large electron-phonon 
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interaction and also the oscillator constant should not be too large. Finally, if we fix g 
and K, the situation for a is more subtle. We see that it should not be too large and 
also it should not be too small. In particular, we are not allowed to set a = and thus 
we cannot conclude anything concerning the Holstein model with a fiat branch. 

A surprising point is the absence of the mass of the phonons in the derived conditions. 
It would not have been the case if we would have considered the sum rule ( ^.141 ). Indeed, 



in this case, we would have transferred inequality ( [4.5| ) onto the two-point correlation 



function for the Sx variables. For this, we need to compute the double commutator 
[s_fc, [-ffA, Sfc]] = — that now depends on m. With this method, we obtain the inequality 



^(-.o3.o)a > (^^)A-(^/d^4v 4am(2d + 2E.cosA;. 



coth 4 / — ; ^= -— . (5.6) 

^' 16ma(2t/ + 2E.cosA;i) ^ ' 

Unfortunately we are unable to use it since we do not know how to derive a lower bound 
for (s^)a. However, this discussion suggests that a lower bound on (s^.)a would yield 
existence of long-range order for another region in the space of parameters which would 
now depend on m. 

6. The hexagonal lattice 

Using the same technique, we can prove a similar result for a two-dimensional hexag- 
onal lattice at zero temperature. In this case however, we do not need to impose an 
external magnetic fiux on the system. We can thus choose real positive hopping am- 
plitudes txy = 1 for \x — y\ = 1 and in direct space the Hamiltonian is given again by 



the sum of (2^), (2^), (|2.5|) . This Hamiltonian is still invariant under a particle-hole 
transformation and we have {nx)A = \-, {sx)a = 0. The two long-range orders we are 
interested in are the ones corresponding to a "chessboard" configuration for the s^ or 
the fix- 

The hexagonal lattice is composed of two sublattices A and B. In figure |l|, we have 
represented the sites of A as black dots and we have defined the vectors ui, U2, u^ and 
Ui, U2 and f/3. It is natural to define the analog of the Fourier transform of a local 
observable Ox using a basis of eigenfunctions of the Hamiltonian ^^ ^^ ^xyC^Cy for a 
single particle hopping on the hexagonal lattice. The eigenfunctions and eigenvalues 
are parametrized by a wave vector k in the first Brillouin zone of the triangular lattice 
defined by A (one has to notice that this triangular lattice has |A|/2 sites and thus we 
have also |A|/2 fc's in the first Brillouin zone) and by another quantum number r = ±1. 
The eigenenergies are given by e^^r = 'tEu with 

Ek = V3 + 2(cos(A; ■ Ui) + cos(A; ■ U2) + cos(fc ■ U^)), (6.1) 

and the eigenfunctions are 

^kAx) = -^=, xe A 



g-ifc-x 






7 

-ik-x 



^kA^) = '^ , X G 5 (6.2) 
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Figure 1. The hexagonal lattice. 



where 



A*=5: 



■^ pife-Mj 



i=i 



Eh 



They satisfy the orthogonality relation 

We define the Fourier transform on the hexagonal lattice as 



Ofc,, = ^^fe,,(x)0^. 



xeA 



(6.3) 



(6.4) 



(6.5) 



The long-range orders we are interested in are the one in the two-point correlation 
functions {s^k,TSk,T)A and {{n^k,T — l/2)(^A:,r — 1/2))a at A; = and r = — 1. The reason 
is that at A; = and r = —1, we have 



So,-i 



mE^--^^E 



xGA 



|A|f 



(6.6) 



GB 



since Aq = 1. Long-range order in this last operator will indicate the presence of a 
chessboard configuration. We will prove that there exist ei and €2 independent of A such 
that 

TT-\{so,-lSo,-l)A = -r^{C}2s^-^S^)C^Sy-'^Sy))A ^ Ci > (6.7) 



xeA 



xeB y<=A 



yeB 
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and 

■r^(po -iPo -i)a = 

-i^((^(n. - 1/2) - 5^K - 1/2))($^K - 1/2) - $^K - 1/2)))a ^ e^ > 

' ' xSA zS-B j/SA J/6-B 

(6.8) 

for large A and zero temperature. The global strategy is exactly the same as for the 
square lattice in two dimensions. We prove the coexistence of the two long-range or- 
ders, find upper bounds for the Duhamel two-point functions, then for the two-point 
correlation functions and finally we use the sum rule 

E E (P-'^,rPk,r)A = E(K - 1/2)^)a = ^ (6.9) 

k T=±l xeA 



to conclude the proof. To prove the first equality in (|6.9| ), one uses the fact that Xk^-k = 
Xk^l = 1. Let us now indicate the main steps of the proof insisting on the points that 
differ from the square lattice case. 

Coexistence of the two long-range orders 

For the coexistence of the two long-range orders and the relation between the Duhamel 
two-point functions, we have the commutation relation 

[HA,p.] = iS^ + iUiny-l/2) (6.10) 

where 

3 

i=\ 
3 

S^ = {KWla)s^\^a^s^_uv x e B. (6.11) 

1=1 
In Fourier transform, we have 

Sk,r=[^ + M^ + rEk)]Sk,r (6.12) 

and for the point A; = 0, r = — 1 in which we are specially interested, Sq-i = ksq-i. 
With these facts, it is easy to prove the coexistence of the two long-range orders ( |6.7|) and 
( |6.8|) . Then, by (|3.6|) , we obtain the relation between the Duhamel two-point functions 



[k + 4a(3 + tEj,)] J {s^k,T, Sfc,r)A - UJ {p^k,T, Pk,T 



^ ^K + 4ai3 + rEk)]. (6.13) 
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Infrared bounds 



To derive an infrared bound for {S-k,T, Sk,T)A we proceed as in eqs. (|4.2|)-P75|). For 
the hexagonal lattice one can bring ( |4.1| ) into a reflection positive form when the flux 
through each plaquette is zero. The analog of (|4.3|) is 



x£A i=l y£A j=l xGA i=l 

(6.14) 
We now choose hx^x+u for a; G A to be (for each k in the first Brillouin zone and each r) 

hx,x+u, = -^ e-''^-^ (1 + T\k e-'^-^^ ). (6.15) 



Inserting this form into (|6.14|) , we get 



We can transfer this information onto an inequality on the Duhamel two-point functions 
for the electrons using (|6.13|) : 



gy{p-k,T,pk,T)A ^ 



K + 4a(3 + rEk) / /t + 4a(3 + rEk) 



Aaf3{3 + rEk) V ^ 



leading to 



K + 24a 
^ , 6.17 

v/a/3(3 + rEk) 



g^a/3{3 + rEk) p 



Inequalities on two-point correlation functions 

As for the square lattice, the upper bound for the two-point correlation function in 
the electronic density is given by 



(p_fe,rPfe,r)A ^ -^/Bk,rCk,rCOth\ ^'^ (6.19) 

2 y ^i^k,T 

where B^^t is the bound given in ( |6.18| ) and Ck,T is the upper bound on the double 
commutator 



(Kfe,[i7,nfc]])A = i^ 5^ t,.,(2-re-''=-(^-^)Afe-re''=-(^-^)A_,)(4c,)A ^ 12 = Ck,r. 

(6.20) 



lAI ^ 
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Finally, in the zero temperature limit (3 ^ oo and for large volumes, we get 



Evaluating the integrals in ( |6.20| ) leads to the final condition on the hexagonal lattice 
for zero temperature 

K + 24a 

^ < 0.18. (6.22) 

g^/a 

The qualitative behavior in the parameters g, a and n is the same as for the square 
lattice. 

7. Conclusion 

In the present work we have provided a rigorous proof of the stability of the Peierls 
instability at half-filling when quantum and thermal fiuctuations are taken into account. 
Our results are valid for dimensions greater or equal to two for a modified Holstein 
model. The regime investigated here is limited to strong coupling but is independent of 
the strength of the quantum fiuctuations. Of course, we expect that for weak coupling 
the mass m of the oscillators should play a role but our method does not give any 
information about this situation. An important drawback of infrared bounds used here 
is that we cannot treat the one dimensional case. 

One would also like to treat the standard Holstein model where the dispersion relation 
is fiat. However we do not know how to obtain the appropriate infrared bound because it 
is not apparent how to couple the phonon field s^ to the h^y variables. Another strategy 
would be to use so-called chessboard estimates to perform a sort of Peierls argument. 
We have carried out this program for a simplified version of the Holstein model where 
the Einstein oscillators at each site x G A are replaced by two levels systems modeled 

(3) 

with Pauli matrices. More precisely, in the Hamiltonian, Sx is replaced by ai and 
p^/2m by eax (e ~ 1/m). For such a model, it is found that in two dimensions the 
Peierls instability occurs for all " electron-phonon" coupling for weak enough quantum 
fiuctuations e [p^ . 

Another important issue is the effect of electron spin and electron-electron Coulomb 
interaction. Most of the analysis carried out here is still valid if we consider spin and 
add to the Hamiltonian a Hubbard term U J^xeAi^^l " V2)('^ii — 1/2)- The coexistence 
theorem is still valid when applied for Sk and Pk = jh X]xeA(^^T + ^^^1 ~^)- ^l^o since the 
Hubbard term is reflection positive at half-flUing we may deduce again infrared bounds 
for (s_fc, Sfc)A and {p^ki Pk)k- The extra information that needs to be controlled is a lower 
bound for 

^{p-kpk)k = 5^(Kt + nxi - 1)^)a = 2 '^{nx^nxi)A. (7.1) 

k xeA xgA 

In other words, one needs to prove that there is a flnite density of doubly occupied sites 
in order to obtain a Peierls instability. When ?7 — > oo, {nx-\nxi)A ^ so that we do not 
expect to have a Peierls instability and therefore it should be present only for moderate 
Coulomb interaction. Let us notice that if (|7.1| ) is of the order of 7|A| for some 7 > 0, 
then the analysis presented in this work automatically implies that at half-fllling there 
is a CDW period two state. 
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